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ABSTRACT 

On 

We calculate the 0(a ew mf, b ,/m^ v ) and 0(a ew mj, b ../m^ v ) supersymmetric elec- 
troweak corrections to the process H + — ► W + h in the Minimal Supersymmetric 



Model. These corrections arise from the virtual effects of the third family (top and 



Q," bottom) quaurks and squarks (top-squark and bottom-squark) . We find that for 

m H + > 200GeV at low tan/3(< 3), the corrections can increase the tree-level decay 



widths and the branching radios more than 20% and 40%, respectively. 
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1. Introduction 

The minimal supersymmeytic standard model(MSSM) takes the minimal Higgs 
structure of two doublesfl], which predicts the existence of three neutral and two 
charged Higgs bosons h,H,A, and H ± . When the Higgs boson of the Standard 
Model(SM) has a mass below 130-140 Gev and the h boson of the MSSM is in the de- 
coupling limit (which means that H + is too heavy anyway to be possibly produced), 
the lightest neutral Higgs boson may be difficult to distinguish from the neutral Higgs 
boson of the standard model(SM). But charged Higgs bosons carry a distinctive sig- 
nature of the Higgs sector in the MSSM. Therefore, the search for charged Higgs 
bosons is very important for probing the Higgs sector of the MSSM and, therefore, 
will be one of the prime objectives of the CERN Large Hadron Collider (LHC). In 
fact, in much of the parameter space preferred by MSSM, namely M#± > Mw± 
and 1 < tan/5 < m t /m6 [2,3], the LHC will provide the greatest opportunity for the 
dicovery of H ± particles. Previous studies [4,5] have shown that the dominante de- 
cay modes of the charged Higgs boson are H + — > t + b with m#± > m t + m^ and 
H + — ► t + v. However, recent analyses[6,7] indicate that the altermative decay chan- 
nel H + — ► W + h° could be very important. In fact, its branching ratio can be rather 
large, competing with the top-bottom decay mode and overwhelming the tau-neutrino 
one for M#± > m t at low tan/3[7]. Thus an more accurate calculation of the decay 
mechanisms is also necessary to provide a solid basis for experimental analysis of ob- 
serving H + — ► W + h° at the LHC. In Ref.[6], R. Santos et al. calculated the top quark 
loops corrections to the decay width of the process H + — ► W + h° in the framework 
of the two-Higgs doublet model for some reasonable choice of the free parameters 
and found that such corrections can be as large as 40%. In this letter, we present 
the calculation of the 0(a ew ml^/rri^ v ) and 0{a ew rnt h Jrn\ v ) supersymmetric(SUSY) 
electroweak(EW) corrections to this process in the MSSM. These corrections arise 
from the virtual effects of the third family (top and bottom) quaurks and squarks 
(top-squark and bottom-squark) . We will give attention manily to these corrections 
in parameter range tan/5 < 3 and m#+ > 200 GeV, where the H + — > W + h° decay 



rate is significant. 

2. Calculations 

Feynman diagrams contributing to supersymmetric electroweak corrections to 
H + — > W + h are shown in Fig.l and 2. We carried out the calculation in the t'Hooft- 
Feynman gauge and used dimensional reduction, which preserves supersymmetry, for 
regularization of the ultraviolet divergences in the virtual loop corrections using the 
on-mass-shell renormalization scheme [8], in which the fine-structure constant a ew and 
physical masses are chosen to be the renormalized parameters, and finite parts of the 
counterterms are fixed by the renormalization conditions. The coupling constant g is 
related to the input parameters e, mw, and mz via g 2 = e 2 /s 2 u and s^ = l — m^/m 2 ?. 
As far as the parameters (3 and a, for the MSSM we are considering, they have to be 
renormalized, too. In the MSSM they are not independent. Nevertheless, we follow 
the approach of Mendez and Pomarol [9] in which they consider them as independent 
renormalized parameters and fixed the corresponding renormalization constants by 
a renormalization condition that the on-mass-shell H + lui and hll couplings keep the 
forms of Eq. (3) of Ref . [9] to all order of perturbation theory. 

The relevant renormalization constants are defined as 

m 2 W0 = m 2 w + 5m 2 w , ra| = w| + Sm 2 z , 

tan/3o = (1 + SZp) tan/?, 

sin a = (1 + SZ a ) sin a, 

W^ = (1 + 5Z w f/ 2 W^ + iZ)H w± d»H*, 

H± = (1 + 6Z H ±) 1 ' 2 H ± , 

H = (l + 5Z H ) 1 / 2 H + Z 1 Jlh, 

h Q = {l + 5Z h ) l ' 2 h + Z l ^H. (1) 

Taking into account the 0{a ew rn 2 n b Jrri 2 v ) and 0{a ew rnt h Jrn^ ¥ ) SUSY EW cor- 
rections, the renormalized amplitude for H + — ► W + h can be written as 

M ren = M + 5M, (2) 



where Mq is the tree-level amplitude arising from Fig. 1(a), which is given by 

Mo = -|cos03-a)p%(£;). (3) 

Here p is the momentum of the incoming charged higges Boson, and k is the momen- 
tum of the outgoing W + Boson. SM represents the SUSY EW corrections, which is 
given by 



l JL C os(P-a)p^(k)fe 

^ 9 



5M = — — cos(/3 — a)p^e^{k)[ tan(/3 — a) (cos (5 sin j35Zp — tdna5Z a ) 



+^(5Z H+ + 5Z h + 5Z W ) - tan(/3 - a)5Z]H] 
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+igm w sin(/3 - a)Z l J^ + i ^ A^e^k), (4) 



with 



5g 5e 1 5m 2 z 1 5m 2 z — Sm^ 
g e 2m| 2 m|- m 2 w 

Sg l tfff 1 »% 7 i/2 

^ " "7 2 <T " 2^ + " ^n7 W ' 

<$# lorajv 1 1/2 . 2 

oZ Q = h - — 5 -dZ h - cot az^ - sin ^043. (5) 
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Here Aj(j = 1-4) are the vertex form factors coming from Fig.l(6)-(e). The 5e/e 

H{b)i 



appearing in Eq.(5) does not contain the 0(a ew m 2 ( b \/m 2 v ) corrections and needs not 



be considered in our calculations. 

Calculating the self-energy diagrams in Fig. 2, we can get the explicit expressions 
of all the renormalization constants as following: 
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-4h b h t m b m t PiifoiB' H+bt + ]T (^(^C^'i + ^'if^o 
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*Z fc = ^b{2^a? 1 (Bf" + 2m?*?" + m 2 <<) + 2^ 1 (B« 



_2 |r)^2^2 i riliit i o..,2r>'/'» .,:'! -n'liit\ i o/,2, % 2 / r; /<'<'< 



+(^^e? yi ) 2 B^]}, 

W = - ^^ r [(M»«0ii + h b m b p 12 )((ml - m 2 )(B° bt - B ^*) - m^sf" 6 *) 

+ E tt^i^^fyi + M>WK - m|)(^^ - Bf~ bii % 
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1fi 2 3 /T 12 2 , »K(i + ^) + ^ 2 -|(i + < A )], 



(6) 



with 



■. s , A r 1 mfyiy — 1) + m 2 (l — y) + mly„ , s 

£f = "1 H— rr " / c%/ n hiK^ T }' 7 

n + 1 Jo fi z 

B i jk = (-1) n / d — y — \_y) (8) 

Here h b = gm b / \/2m w cos/? and /i t = gm t / \/2m w sin /3 are the Yukawa couplings 
from the bottom and top quarks, respectively. The notations #*-, 6*^-, «„, /3y, </?y and 
Ofj k used in the above expressions are defined in Appendix A. 

Calculating the diagrams in Fig.l(6)-(e), we can get the explicit expressions of the 
vertex form factors Aj(j = 1-4) as following: 

Ai = ~^ sma{h b sinP(Bf bt - Bf tb ) + [-h t m b m t cos/3(C + 2C 2 ) 
-2h b sin (3C 00 + hm 2 b sin p(C - 2C 2 ) - h b m 2 H+ sin/3(C 12 + AC 22 ) 
+h b m 2 h sin/5(Ci2 + 2C22) + h b m 2 w sin/3(Ci - 3C± 2 



-4C 22 )] (m 2 w , m 2 H+ ,m 2 h , m 2 b , m 2 , ml)}, 
A 2 = TF% cos a{h t cos(3(-Bg +bt + B* +ht ) + [h b m b m t sin/3(C + 2C 2 ) 

l07T z 



+2/it cos/3C o - to* cos/?(C - 2C 2 ) + ^m^+ cos/3(Ci 2 + AG 



22 J 



+h t m 2 h cosl3(Ci2 + 2C 22 ) - /i t m^ cos/3(Ci - 3C 12 
-4C 22 )] (m^, m 2 H+ ,m 2 h , m 2 , m 2 , m 2 )}, 

As = g To 2 5Z 51 ^ii^u'i^jj') 9ki9lk'®i'j'i( h b®j>k'i + h&j'k'i) 

xC 2 (m 2 v ,m 2 H +,m 2 h ,ml i , m \ > m | ) ' 

A4 = ~ Q To 2 5Z 51 ^kl^kk'^ii'i^jf) ®j>i'l(hb®k'j'l + h t® k ij>l) 

xC 2 {m 2 Wl m 2 H+ ,m 2 h , m 2 ., mf fc , m|), (9) 

where Cj and C^ are the three-point Feynman integrals [10]. 
The corresponding amplitude squared is 

\M ren \ 2 = \M \ 2 + 2Re(5MM t ). (10) 

The renormalized decay width is given by 

T(H + -»■ W + /i) = r (# + -»• W + h) + 5r(iJ + -► W + h) 

= 77, —\/( m l - m H+ ~ m w) 2 ~ 4m^+m^|M ren | 2 . (11) 

16nm 6 H+ v 

The tree-level branching ratio of if + — ► W + h decay is 

b (h + - py+/i) = r (#+ - w + / i )/r (# + ), (12) 

where the tree- level total decay width To(H + ) of the charged Higgs boson is approx- 
imated by 

r (# + ) = T (H + -► W + h) + T (H + -► tb) + Y (H + -> cs) + T (iJ + -► ri/). (13) 

While calculating the one-loop branching ratio B{H + — ► VT + /i), the QCD correc- 
tions [11] to the i/ + — ► £6 width, <5r(i/ + — ► £&), are included into the total width 
T(H + ), but its leading EW corrections[12] are neglected since they are much smaller 
than the QCD corrections: 

B(H + -► W + h) = T(H + -► W + /i)/r(#+) (14) 
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with 

T(H + ) = T {H + ) + 5T(H+ -► W + h) + 5T{H + -+ £6) (15) 

3. Numerical results and conclusion 

We now present some numerical results for the SUSY EW corrections to the 
decay H + — ► W + h. The SM input parameters in our calculations were taken to 
be a ew {m z ) = 1/128.8, m w = 80.375GeV and m z = 91.1867GeV[13], and m t = 
175.6GeV and m^ = 4.7GeV. Other parameters are determined as follows 

(i) The one-loop relations[14] between the Higgs boson masses M hH ^ H T and the 
parameters a and 8 in the MSSM were used, and m#+ and 8 were chosen as the two 
independent input parameters. As explained in Re [9], there is a small inconsistency 
in doing so since the parameters a and 8 of Ref[14] are not the ones defined by 
the conditions given by Eq.(Eq(3)) of Ref[9]. Neverthless, this difference would only 
induce a higher order change [9]. 

(ii) For the parameters mj~ - ^ and A t ^ in squark mass matrices 

/ Mh m„M LR 
Ml=\ LL " | (16) 

\ m q M RL M 2 RR 

m Il = m Q + m l + m % cos 2f3{lf - e q sin 2 6 W ), 

O O O /~t O s\ 

M RR = rrijj p + m q + m z cos 2pe g sin ^^, 

I A t -u cot 8 (q = i) \ 

M LR = M RL = „ , (17) 

y A 6 — /itan/9 (q = b) J 

to simplify the calculation we assumed Mq = M^ = M^ and A t = A^, and we used 
Ma, A t and /i as input parameters except the numerical calculations as shown in 
Fig.5(a)and (b), where we took m^, m~ b , A t — A\ and /x as the input parameters. 

Figure 3 shows the tree-level partial width are relatively large for low values of 
tan/3(= 1.5,2) when m H+ > 200GeV. 



with 



In Figs. 4 (a) and (b) we present the SUSY EW corrections to the tree-level decay 
widths and the braching radios as functions of ra#+ for different values of tan/3, 
respectively, assuming Mq = M- t = M b = 300GeV, A t = A b = 300GeV, and 
\x = — lOOGeV. Figure 4(a) shows that the relative corrections increase the par- 
tial width significantly at low tan/3, which can exceed 40% for tan/3 = 1.5 and 
30% for tan/3 = 2, respectively. Fig. 4(b) gives the tree-level branching ratios and 
the branching ratios after including (^(ae^m^/m^) and O^aewmt^/m^y) SUSY 
EW corrections for tan/3 = 1.5 and 3, respectively. From Fig.4(b), we see that the 
branching ratios are enhanced by the corrections, especially for m#+ > 300 GeV, 
they can be increased by 50%. But the corrections to the branching ratios decrease 
with an decrease of mn+ , especially for m#+ bellow 200 GeV, they become negligibly 
small. 

In Fig.5(a) and (b) we assumed Mq = M h A t = A b = 500GeV, \x = lOOGeV, 
m^ = 170GeV and m b _ = 200GeV. Fig. 5 (a) shows the relative corrections to the tree- 
level widths as a function of m#+ for tan/3 = 1.5, 2, 6, 10 and 30. Comparing Fig. 5(a) 
with Fig. 4(a), we see that the relative corrections are enhanced by the relatively light 
squark masses in the case of heavy charged Higgs masses at low or high tan/3. Fig. 5(a) 
also shows at low tan/3 the corrections always increase the decay width , which can 
exceed 50%, while at high tan/3 the corrections decrease the decay width significantly. 
There are dips at m#+ = m^ + m b = 370GeV due to the threshold effects. As shown 
in Fig.4(b), Fig.5(b) shows the branching ratios are enhanced by the the SUSY EW 
corrections for tan/3 = 1.5 and 3, especially when m#+ > 300GeV. 

In conclusion, we have calculated the O^e^m^/m^/) and O^a^rn^^lrn^) SUSY 
EW corrections to the process H + — > W + h in the MSSM. In general, these corrections 
increase the decay widths and the branching radios and are sensitive to both of tan/3 
and the mass of the charged Higgs boson. For m#+ > 200GeV at low tan/3(< 3), the 
corrections can increase the decay widths and the branching radios more than 20% 
and 40%, respectively. 
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Appendix A 

We present some notations used in this paper here. We introduce an angle <p = (3 — a, 
and for each angle a, (3, <p, 9 l or 9 b , we define 



O; 



9 h 



cos a sin a 

— sin a cos a 

cos 9 l sin 9 l 

— sin 9 l cos 9 l 



i Pij 






cos (3 sin f3 

— sin j3 cos (3 

cos 9 h sin 9 b 

— sin 9 b cos 9 h 



<Pij 



cos v? sin (p 



sin </? cos (/? 



We define six matrix e i j kl (i =1 — 6) for the couplings between squarks and Higgses: 



B 1 



el 



i.)2 



b 2 



B 2 



b 5 

u ijl 



B 5 



B 6 

u iii 



2m t cos ct A t cos a + /x sin a 

A t cos a: + fi sin a 2m 4 cos a 

2m t sin a A t sin a — fi cos a 

A t sin a — fi cos a 2m t sin a 

2rrib sin a Ab sin a + fi cos a 

ylb sin a + /i cos a 2rrib sin a 

Irrtb cos o; A5 cos a — fi sin a 

^4;, cos a — fi sin a 2m;, cos a 

mb sin (3 

^ ^ sin (3 + /i cos /3 m t sin /3 

— nib cos P 

— A fe cos /3 + fi sin /3 

m t cos/3 A t cos /? + fi sin /3 
m fe cos (3 
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o 6 



m t sin /3 A t sin (5 — \i cos (3 
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Figure Captions 

FIG. 1 Feynman diagrams contributing to supersymmetric electroweak corrections 
to H + — ► W + h: (a) is the tree-level diagram; (b) — (e) are the one-loop diagrams. 

FIG. 2 Feynman diagrams contributing to the renormalization constants. 

FIG. 3 The tree- level H + -»■ W + h decay widths. 

FIG. 4 (a): The SUSY EW relative corrections to the H+ -> W+/i decay width; 
(b): The branching ratio of the H + — > lU + /i decay assuming Mq = M^ = M^ = 
300GeV, A t = At> = 300GeV, and \x = — lOOGeV. Logarithmic coordinate is 
used in (b). 

FIG. 5 (a): The SUSY EW relative corrections to the H + — > W + h decay width; 
(6): The branching ratio of the H + — ► IU + /i decay, assuming Mq = Mf, A t = 
A 6 = 500GeV, /i = lOOGeV, ra ?1 = 170GeV and m~ bl = 200GeV. Logarithmic 
coordinate is used in (b). 



13 



H+- 



(a) 



H+- 



>W+ 



h 




ty 
< 

(d) 



W 



H+- 



h 



V 
< 

(e) 



W + 



h 



Fig.l 





H+ 




H+ 



t 
(c) 



---4 f--- 

t 

(d) 



H^ 




---4 



/ w+ 



t 



t(b) 




h V J h ( H ) 
t(b) 

(a) 



1(b) 



h \ / h (H) 



i(b) 
(h) 



Kb) 

/ \ 

[ I 

h h(H) 

(0 



Fig.2 



Fig. 3 



0.05 



0.04 - 



> 
CD 

O 



0.03 - 



U 0.02 - 



0.01 - 



0.00 




150 



200 



250 



300 350 

m H+ (GeV) 



400 



450 



500 



Fig. 3 (a) 



CO 




150 



200 



250 



300 350 

m H+ (GeV) 



400 



450 



500 



Fig.3 (b) 



10 U r 



05 
DC 



O 

c 

05 



10- 7 



00 10 2 



B (tanp = 1.5) 
B(tanp = 1.5) 
B (tanp=3) 
B(tanp=3) 




150 



200 



250 



300 350 

m H+ (GeV) 



400 



450 



500 



Fig.4(a) 



CO 



o 20 - 




■40 - 



200 



250 



300 350 

m H+ (GeV) 



500 



Fig.4 (b) 



10 U r 



rt 10 1 
DC 



o 

c 

,-2 



10" r 



B (tanp = 1.5) 
B(tanp = 1.5) 
B (tanp=3) 
B(tanp=3) 




150 



200 



250 



300 350 

m H+ (GeV) 



400 



450 



500 



